Momentum distribution and Compton profile by the ab initio GW approximation 
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We present two possible approaches to calculate the momentum distribution n{p) and the Comp- 
ton profile within the framework of the ab initio GW approximation on the self-energy. The ap- 
proaches are based on integration of the Green's function along either the real or the imaginary axes. 
Examples will be presented on the jellium model and on real bulk sodium. Advantages and draw- 
backs of both methods are discussed in comparison with accurate quantum Monte Carlo (QMC) 
calculations and x-ray Compton scattering experiments. We illustrate the effect of many-body cor- 
relations and disentangle them from band-structure and anisotropy effects by a comparison with 
density functional theory in the local density approximation. Our results suggest the use of Go Wo 
momentum distributions as reference for future experiments and theory developments. 



PACS numbers: 71.10.-w, 78.70.Ck, 71.20.Dg, 02.70.Ss 

INTRODUCTION 

The momentum distribution n(p), defined as the prob- 
ability to observe a particle with momentum p, is one 
of the basic quantities in quantum statistical mechanics. 
For fermions, it is a quantity showing a direct evidence 
of the Pauli principle (see Fig. [T]) . For an ideal Fermi gas 
of free non-interacting electrons, the momentum distri- 
bution is the Fermi-Dirac distribution. In equilibrium at 
zero temperature, it is the step function, i.e., it is 1 for 
p below the Fermi momentum pp and above, with a 
discontinuity equal to one at the Fermi sphere surface. 

When passing from the ideal Fermi gas to a Fermi liq- 
uid of interacting fermionsji the momentum distribution 
departs from the perfect step function. Correlations in- 
duce a modification of the distribution with a spill out 
from lowest to highest momenta, so that the probability 
to observe an electron at a momentum p larger than the 
Fermi momentum pp becomes finite even at zero tem- 
perature. However, n{p) still retains the discontinuity at 
p = Pp, having only its magnitude reduced from 1 . 

For the homogeneous electron gas (HEG) or jellium, 
one of the most fundamental models to study electronic 
correlations, the discontinuity is reduced to the quasipar- 
ticle renormalization factor Zpp calculated at the Fermi 
momentum.^ In the high-density limit (ts = 0), which 
is dominated by the kinetic energy, Zpp approaches the 
uncorrelated value 1. The renormalization factor Zp^, 
is expected to reduce with decreasing density as corre- 
lations build up (see Fig. [T]). The discontinuity is still 
retained at finite densities and vanishes for oo. In 

particular, for so-called strongly correlated systems, the 



discontinuity is strongly suppressed. The modification 
of the momentum distribution and the reduction of the 
discontinuity is mainly a correlation effect, unaffected by 
other, Hartree or exchange, many-body effects. Thus, 
the momentum distribution and its discontinuity, unlike 
other observables such as the band dispersion or the gap, 
provides a unique and unambiguous quantification of the 
level of correlation in a system/2."J-3 Therefore, experi- 
mental measurements of n{jp) and its discontinuity, are 
of fundamental importance to test and verify many-body 
theories. The level of accuracy in describing correlations 
by a given theoretical approach can be directly quantified 
by a comparison of calculated n(p) and Z with experi- 
mental results, if they were available. 

In this work, we focus on the calculation of the mo- 
mentum distribution and related quantities within the 
framework of ab initio many-body theory in the GW ap- 
proximation. We have studied two possible approaches: 
the momentum distribution can be calculated by inte- 
gration of expressions containing the Green's function or 
the self-energy evaluated on either the real or the imag- 
inary uj axis. Analytically, both integrals provide the 
same result. Numerical convergence problems with re- 
spect to the integration sampling favour the real-axis in- 
tegration for evaluating the momentum distribution far 
away from the Fermi momentum, while the imaginary 
axis is more accurate near the discontinuity. Taking bulk 
sodium and jellium as examples, our results show that 
the G{)W{) momentum distributions are in good agree- 
ment with x-ray Compton scattering experiments^'* and 
also with quantum Monte Carlo results^, provided that 
the appropriate methodology, as analyzed in this article. 



is used in their calculation. As standard GqWq are ap- 
plicable to describe a broad range of realistic systems, 
results on the momentum distribution and related quan- 
tities, e.g. the quasiparticle renormalization factor or the 
Compton profile, can provide accurate reference values 
and stimulate new experiments. In turn, since experi- 
mental measurements of the Compton profile derivative 
discontinuity provide an unambiguous quantification of 
the level of correlation, the limit of validity of the GW 
approximation can be checked, e.g. in strongly-correlated 
systems, stimulating further theoretical progress. 

The momentum distribution n{p) is also the Fourier 
transform of the density matrix n{r, r') which is the fun- 
damental degree of freedom of density-matrix functional 
theory (DMFT)i^"ii. Since DMFT is an in principle ex- 
act theory to calculate the density matrix, it can give 
access not only to ground-state observables such as Kohn- 
Sham DFT, but also to some spectral information such as 
the quasiparticle renormalization factor Z at the Fermi 
surface, which can be exctracted from the discontinuity of 
n{p), for example by switching on and off the pseudopo- 
tential or the electronic correlations. Momentum distri- 
butions obtained within GW, as presented here, can also 
represent a good starting point for improving exchange- 
correlation approximations within DMFT. 

In previous worksfi^"— Kubo has addressed the mo- 
mentum distribution in solids by a GqWq real axis inte- 
gration method using the Hamada plasmon-pole model^i 
which has been criticized^^ for not providing an accurate 
description of the imaginary part of the self-energy, hence 
of the momentum distribution. Here we calculate and use 
the full frequency dependence of the dielectric function. 

The paper is organized as follows: Section|T]is an intro- 
duction to the momentum distribution. In SecHIl we dis- 
cuss the relation with the experiment. In Sec lIIIi we dis- 
cuss periodic crystalline band-structure and anisotropy 
effects on n(p). In Sec. IIVI we consider purely many- 
body effects and present details of the two GW method- 
ologies used in the calculation. Finally, we draw conclu- 
sions and an outlook. Unless explicitly specified, we use 
atomic units (a.u.). 



I. THE MOMENTUM DISTRIBUTION 

The spin averaged momentum distribution is defined 

as 

(J cr 

where ^' is the ground-state wave function of the system 
and fipa = a|,g.apcr is the product of creation/annihilation 
operators for electrons at momentum p and spin a. By 
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FIG. 1. The GqWo momentum distribution of the jellium 
model at decreasing densities, as an illustration of increasing 
correlation effect. Gray short-dashed line: ideal Fermi gas; 
black long-dashed line: rs = 1; light-green dot-dashed line: 
rs = 2; blue dotted line: Vs — 4; red solid line: = 5; 
dark-green dashed line: = 10; 

replacing, in terms of the field operators, ■i/;J.(r), V'o-Cr), 
^(P) = ^E-^ / drdr'e-'P("-^')(vI.|V't(r')V'.(r)|*) 

= drdr'e"'P('"-''')n(r,r'), (1) 

the momentum distribution is expressed as the Fourier 
transform of the reduced single particle density matrix 
n{r,r'). One can see that density-matrix functional the- 
ory is particularly suitable to calculate the momentum 
distribution n{p) and its discontinuity C, related to the 
quasiparticle renormalization factor Z at the Fermi sur- 
face. 

The momentum distribution is the step function only 
in the ideal case of a system of non-interacting free 
fcrmions in equilibrium at zero temperature. The depar- 
ture from that ideal shape is due to many factors that 
must be taken into account, namely: 

A) finite temperature effects; 

B) band-structure effects; 

C) anisotropy of the Fermi surface; 

D) electron-phonon interaction effects; 

E) electron-electron correlation many-body effects. 

In the next sections we estimate and discuss the weight 
of each effect compared to the others. 

In the paper we use the convention to define the discon- 
tinuity of n{p) as C, to be distinguished from the quasi- 
particle renormalization factor Z at the Fermi surface, 
and affected by all other, e.g. band-structure etc., effects. 
It should be noted that ( = Zp^ only holds in the case 
of jellium. 
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A. Non-interacting free fermions and 

finite-temperature effects 

Consider an unpolarized non-interacting and free elec- 
tron gas of given density n = S/Anr^ and Fermi mo- 
mentum pf = {dn/Ay^^/rs- The electronic structure 
is simply provided by perfectly parabolic energy levels, 
Sp = and plane- wave orbitals, (/)p(r) — e*P ''/v^- 

Due to their fermionic nature, electrons obey the Pauli 
principle, and at equilibrium single-electron states are 
filled according to the Fermi-Dirac distribution at tem- 
perature T and chemical potential fi with occupation 
numbers n{p) = l/(e'-^ /2-m)/t _|_ which in this case 
coincides with the momentum distribution. In the limit 
T — > 0, n{p) reduces to the step function, n{p) = 
6{pF ~ p), with a well-defined jump at pp of magnitude 
C — n{pp) — n{p'^) — 1 (gray dashed line in Fig. [IJ. 
Finite temperatures introduce a smoothing of the jump 
in the region — p'^\/2m < T. At room temperature, 
\P/PF - 1| < T/Tf where Tp = p|,/2 is the Fermi tem- 
perature. For metals, we typically have T/Tp « 10~^ at 
room temperature, so that temperature effects are fairly 
negligible for our analysis. 

B. Non-interacting electrons in solids: 

Band-structure effects 

In order to describe non-interacting electrons in a pe- 
riodic crystalline solid, we have to use Bloch energies e^]^ 
and wave functions, 

</>.k(r)=^^Sce^(''+°)^ (2) 

G 

where the summation runs over reciprocal lattice vectors 
G, and k and v are the crystal momentum and the band 
index, respectively. 

At T = 0, all states {ly, k} are filled up to the Fermi en- 
ergy ep-, with a step function for the occupation number, 
n^k = S(ep — e^k)- As a consequence, the momentum 
distribution in the independent particle approximation 
becomes 

n(p) = ^ Oiep - 6,k) ^ <5(p - k - G) l^^kp, (3) 

i/k G 

and will deviate from the ideal-gas step function. 

Whereas the Fermi momentum is no longer a well- 
defined concept, a discontinuity of the momentum dis- 
tribution still occurs whenever the Bloch energies cross 
the Fermi surface, defined hy ep. Eventually, more than 
one discontinuity can be present in the momentum dis- 
tribution when the Fermi surface is organized in sev- 
eral branches. The contributions of momentum den- 
sity centered at reciprocal lattice vectors G 7^ are 
called high-momentum ( umklapp ) components of the mo- 
mentum density and can also be seen experimentally in 
non-jellium-like systems.^^ Assuming a partially filled, 



isotropic valence band, the Fermi sphere of which is en- 
tirely contained within the first Brillouin zone (IBZ), 
e.g., the case of Na, the value of the discontinuity at the 
Fermi surface is reduced compared to the free-electron 
gas by the value 

c = \~€^i\A''<^ (4) 

With respect to the electron gas, band-structure effects 
thus reduce the discontinuity to a value provided by the 
coefficients of the plane-wave expansion of the orbitals, 
Eq. . The more the system approaches the ideal Fermi 
gas, the closer the wave functions are to plane waves, and 
the larger is the discontinuity at the Fermi surface. 

Band-structure effects can be calculated by a simple 
single-electron theory. We will discuss in detail in Sec. lIIII 
how to account them by density functional theory and 
will provide an illustration of these effects. 



C. Non-interacting electrons in solids: Anisotropy 
effects 

In a real solid, the momentum distribution defined in 
Eq. ([3]) is in general anisotropic, and the Fermi surface is 
nonspherical. Whereas angle-resolved momentum distri- 
bution measurements^"*"^* may eventually resolve the full 
three-dimensional (3D) n(p), Fermi surface anisotropy 
may prohibit the access to the magnitude C of the n(p) 
discontinuity by a simple powder-averaged measurement. 
We will discuss anisotropy effects also in Sec. IIIII 



D. Electron-phonon interactions 

An electron in a solid can absorb or emit a phonon 
with momentum q and in a mode with a probability 
given by the electron-phonon coupling \^q. The momen- 
tum distribution of electrons turns out to be modified 
as a consequence of such scatterings. In order to pre- 
cisely account for such effects, one should evaluate the 
full electron-phonon scattering matrix and introduce its 
effect into the momentum distribution as a function of 
the temperature. This can be done by constructing an 
electron-phonon self-energy, for example, in the second 
Born approximation. The electron-phonon coupling may 
lead to a further decrease in the discontinuity of the mo- 
mentum distribution. However, since the phonon De- 
bye frequency, ud, is small compared to the Fermi tem- 
perature, changes in the momentum distribution due to 
electron-phonon interactions are expected only within a 
narrow momentum region 5p/pp < ujp)/Tp « 10~^, well 
beyond the current experimental resolution. Similar to 
pure temperature effects on the electronic distribution, 
we will neglect these effects in the following. 
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E. Jellium and Fermi liquid behaviour: 
Electron-electron correlation effects 

In an electron- electron interacting system, collisions 
between electrons will in general reduce the discontinu- 
ity in the momentum distribution. This effect can be 
accounted for only by a theory presenting a good de- 
scription of correlations and in principle exact to calcu- 
late the density matrix n(r,r'), or directly the momen- 
tum distribution n(p). Density- matrix functional the- 
ory can be such a theory, once a good approximation 
for the exchange-correlation functional is found. Den- 
sity functional theory in the local density approximation 
(DFT-LDA) or in the generalized gradient approximation 
(GGA) present a good description of correlations in the 
electronic density n(r) or other ground-state properties, 
but Kohn-Sham DFT is not an exact theory to calculate 
n(r, r') or n{p). 

Precise values of the momentum distribution can be 
obtained from quantum Monte Carlo (QMC) methods. 
These methods are based on the full many-body wave 
functions ^(ri, . . . , tn) and correlations are explicitly in- 
troduced in many-body Jastrow and backflow potentials. 
The reduced single-particle density matrix can be calcu- 
lated by integrating the many-body wave function, 

n{r,r')^N / dra . . . drN **(r, ra, . . . , rN)*(r', ra, . . . , tn) 



Therefore, n(p) is alternatively given by 



duA{r,r',uj). (8) 



(with J dri . . .drj^'^*'^ — 1) and the momentum distri- 
bution is obtained by Fourier transform, Eq. ([1]). Ac- 
curate results on the momentum distribution and on the 
renormalization factor at zero temperature have been ob- 
tained by quantum Monte Carlo calculations ?i2i2£ This 
way is also available for other wave-function-oriented 
many-body theories, such as, for example, quantum 
chemistry approaches. 

Finally, the momentum distribution can also be cal- 
culated in the Green's function approach to many-body 
theory by 

nip) = -^ J drdr'e-P("--') J ^G<(r, r', 

where is the correlation lesser Green's function. For a 
system in equilibrium at T = K, n(p) can be expressed 
as 

n(p) = l f drdr' e-^P^"-""') T duj A{r , r' , u) , (5) 

*^ J J — oo 

where A is the spectral function of single-particle exci- 
tations that can be obtained from the imaginary part of 
the Green's function 



A{uj) = ——^G{ijj) sgn(w — yu), 

TT 



and is normalized as follows 



+ 00 



did A{r, r', w) — S{r — r' 



(6) 



(7) 



Thus, the knowledge of the spectral function A, or of the 
Green's function, can provide the momentum distribu- 
tion upon frequency integration. 

The spectral function of a free, non-interacting gas of 
fermions is diagonal in momentum space, and reduces to 
a delta function in frequency space 

Aip,u:)^S{Lu~py2) 

As a consequence, at zero temperature, where the Fermi- 
Dirac function becomes a step-function, the momentum 
distribution contains a jump at the Fermi momentum, 
pp. The magnitude of the jump is maximal, C — 1: ^-iid 
directly equal to the strength of the delta-function. 

For any approach characterized by an effective, static, 
single-particle Hamitonian, in particular, the Hartree- 
Fock method or the Kohn-Sham scheme to DFT, the 
spectral function is still characterized by exact delta func- 
tions in frequency space. As a consequence, in the homo- 
geneous system, the momentum distribution is still given 
by the step function of the ideal gas with a jump C = 
while in solids the discontinuity is reduced only by other 
effects, e.g., band-structure or anisotropy. 

Correlations will in general lead to a broadening of the 
singic-particle excitation spectra and smooth out the dis- 
continuity in the momentum distribution. However, for a 
normal Fermi liquid, the damping of an excitation at the 
Fermi surface vanishes^ and we approach again a delta- 
peak in the spectral function, but of reduced weight, the 
so-called quasiparticle renormalization factor Z. There- 
fore, we still expect a discontinuity in the momentum 
distribution, but its magnitude is reduced by a factor 
Z < 1. In jellium, where band-structure and anisotropy 
effects are absent, the jump magnitude itself C represents 
a direct measure of the strength of quasiparticle excita- 
tions Zpp at the Fermi surface. From Fig. [I] we expect 
significant changes in the momentum distribution of jel- 
lium at metallic densities compared to free electrons, a di- 
rect signature of electronic correlations in the system. In 
real systems, in particular jellium-like simple metals, the 
measure of the momentum distribution jump can provide 
a direct evidence of Fermi liquid behaviour. For all other 
metals, the measure of once disentangled all other, 
e.g., band-structure, etc., reduction effects, can provide 
an estimate of Z^p and thus a quantification of corre- 
lation effects and possible deviation from Fermi liquid 
behaviour. 



In Sec. lIVl we will enter into the details of a momentum 
distribution calculation in the framework of many-body 
perturbation theory, including correlations by the GW 
approximation. 
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FIG. 2. Compton profile of valence electrons in sodium. Or- 
ange dotted line: experiment; gray short-dashed line: ideal 
Fermi gas; green line: QMC calculation; violet dashed line: 
DFT-LDA calculation; blue squares: GqWq real-axis calcula- 
tion. The ideal Fermi gas result is nearly indistinguishable 
from the DFT-LDA hue. 



II. THE COMPTON PROFILE AND 
EXPERIMENTAL MEASUREMENTS 

An interesting technique that can provide an experi- 
mental indirect access to the momentum distribution is 
Compton scattering of x-rayS)^ a branch of non-resonant 
inelastic x-ray scattering spectroscopy (IXSS) performed 
with hard x rays (with energies 10-100 keV) in the limit 
of large transferred momentum and energy. IXSS mea- 
sures the probability for inelastic x-ray scattering, quan- 
tified by the double-differential scattering cross section 
(Pa/dfldu! in the solid angle element dfl and transferred 
energy interval du). At values of uj much larger than the 
electron binding energy, the experiment is in the regime 
of the so-called impulse approximation^^ , where the dif- 
ferential scattering cross-section is proportional to the 
Compton profile J(q), 

^ « Ji<l): (9) 

where q is the electron momentum vector that points 
in the direction of the scattering vector. The Compton 
profile is directly related to the momentum distribution 
and defined as the projection of n{p) onto the direction 
q of the scattering vector, 

2 f dp 

n J (2ny 

where n = 2 J dp/{2TT)^ "-(p) is the average electron den- 
sity (factor 2 is for spin). 

Substantial simplification is obtained assuming an 
isotropic system n{p) = n{p). In this case, we have 

J^l) = ^ / dppn{p) (10) 



and n{p) is determined by simple differentiation of the 
Compton profile 

2Ti'^n dJiq) 
n{p) = ^ — 

Metals with highly isotropic momentum distribution are 
therefore most suited for measurements of the momen- 
tum distribution via Compton scattering. Notice that 
we are using the convention by which the Compton pro- 
file is normalized to 1, 

j dqJ{q) = l. 

For example, in the case of the ideal Fermi gas the 
momentum distribution n(p) is the step function n(p) — 
9{pp —p) with a discontinuity C = 1 at the Fermi momen- 
tum pp. The associated Compton profile is an inverted 
parabola for q < pp, 

Ap%, 

and vanishes for q > pp (see Fig. ^ . It presents a dis- 
continuity of the first derivative at the Fermi momentum 
q — Pp, which is related to the discontinuity C of the 
momentum distribution n{p). A direct measure of the 
discontinuity of the first derivative of the Compton pro- 
file dJ/dq provides the discontinuity of the momentum 
distribution which we are interested in. 

Measurements of the Cornpton profile of solid sodium 
have been presented in Ref. Il4l . These experiments were 
performed at the beamline ID16 of the European Syn- 
chrotron Radiation Facility. Details of the beamline, the 
spectrometer, sam ple p reparation, and data analysis are 
given in Refs. [13, and [H. To get Compton pro- 
files of valence electrons, the core-electrons contribution 
must be subtracted from the experimental signal. Core- 
electron IXSS spectra were calculated by a quasi-self- 
consistent field (QSCF) approximation^i and the real- 
space-multiple scattering approach FEFFq^^, both found 
in agreement with each other and the measured Comp- 
ton spectra In Fig. [5] we present the resulting experi- 
mental Compton profile of sodium valence electrons to- 
gether with various theoretical results. With respect to 
the Compton profile of an ideal Fermi gas, we observe a 
departure from the perfect inverted parabola and a re- 
duction of the discontinuity in the first derivative at the 
Fermi momentum. In the next sections we will analyse 
in detail the influence of band structure and correlations 
on the Compton profile and the underlying momentum 
distribution within the GW approach. 

III. BAND-STRUCTURE AND ANISOTROPY 
EFFECTS BY DFT-LDA 

Density-functional theory is an in principle exact the- 
ory to calculate the ground-state electron density n(r) = 




FIG. 3. The Fermi surface in sodium, obtained from a Kohn- 
Sham DFT-LDA calculation. 
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(\I'o|V'^(r)V-'(r)|^o) together with the ground-state en- 
ergy. However, neither the off-diagonal elements of the 
density matrix rt(r,r') = (vI'o|V''''(r')V'(r)|\E'o), nor its 
Fourier transform, the momentum distribution, can in 
general be described within DFT, similar to quasiparti- 
cle properties. 

Nevertheless, the Kohn-Sham scheme of DFT can be 
used to evaluate the leading-order band-structure and 
anisotropy effects to the momentum distribution. This 
can be done by replacing Kohn-Sham wave functions and 
energies into Eq. (jS]) where correlation effects are ne- 
glected. Whenever the Kohn-Sham wave functions are 
close to the exact quasiparticle ones, and quasiparticle 
corrections do not significantly modify the Kohn-Sham 
Fermi surface, Eq. ([3|) is expected to provide a good de- 
scription of band-structure and anisotropy effects, as well 
as the position of the discontinuity. This is the case in 
alkali metals and simple semiconductors A breakdown 
of this picture was found in solids which contain shal- 
low d electrons'^^, and may also be expected for systems 
with shallow / electrons. In these cases, the Kohn-Sham 
wave functions of d electrons are in general not close to 
quasiparticle ones. Further, a quasiparticle shift of d lev- 
els may significantly modify the DFT Kohn-Sham Fermi 
surface, and so the position of the discontinuity. 

The system studied here, sodium, is significantly sim- 
pler, and DFT provides a reasonable estimate of band- 
structure and anisotropy effects. We run a standard 
Martin- Trouillers pseudopotential DFT-LDA calculation 
on a plane- waves basis set, with a cutoff of 24 Ryd on 
the kinetic energy, fixing the lattice parameter to the ex- 
perimental value 8.108 a.u. (corresponding to an average 
density of = 3.99). Using a Monkhorst-Pack 8x8x8 
grid of k points to represent the Brillouin zone, the elec- 
tron density is converged at self-consistency. Larger sets 
of k points were used to calculate the band plot, the 
Fermi surface, and the momentum distribution. 

We start with an analysis of the anisotropy effects, 
as provided by DFT. The Fermi surface (Fig. [3]) is very 
close to a perfect sphere with a deviation of 0.2% only. 
This result shows a close similarity between sodium and 
the jellium model. In order to estimate the maximum 



FIG. 4. Anisotropy of the Fermi surface. Sodium first band 
energy vs. k: Blue circles and line: along the F-H direction; 
green squares and line: along the F-F direction; red diamonds 
and line: along the F-N direction. Black line: Fermi gas at 
Ts — 3.9917 energy levels. 



anisotropy of the Fermi surface in sodium, we calculated 
the Kohn-Sham energies on a very fine k-mesh across the 
Fermi level and along the three high symmetry directions, 
F-H, F-F and F-N (Fig.g]). The anisotropy at the Fermi 
surface is evaluated to be Ak = 5.5- 10~^ Bohr"-'^. Such a 
value is well within the accuracy of the IXSS experiment. 
Therefore, anisotropy effects in sodium are currently not 
detectable^, and we do not expect more insights from 
angle resolved measurements on single crystals compared 
to much simpler powder averaged measurements. 

The momentum distribution calculated using Eq. (|3]) 
and Kohn-Sham DFT-LDA wave functions is shown in 
Fig. [5l A very small deviation from the momentum dis- 
tribution of the ideal Fermi gas is observed below and 
above pf- The DFT-LDA discontinuity is evaluated to be 
Cdft = 0.98(1), slightly below the ideal Fermi gas value. 
An analysis of the Kohn-Sham wave functions shows that 
the plane wave coefficients <^^i°k in sodium are 0.99 for 
all k-points around kj?. This value remains almost the 
same for states at the bottom of the first band and up 
to 1 eV above the Fermi level. Therefore, wave functions 
of the first band are very close to plane waves in a wide 
range of fc, consistent with the band-structure (Fig. [6]). 
The first band in sodium is an almost perfect parabola 
which superimposes with the dispersion curve of the ideal 
Fermi gas. A difference between the two curves is appre- 
ciable only above the Fermi level and close to the BZ 
boundary. The DFT band-width is only 0.04 eV larger 
in the ideal Fermi gas than in sodium, a value below the 
standard accuracy of a pseudopotential calculation. 

Beyond the first Brillouin zone, band-structure effects 
introduce deviations from the momentum distribution 
of an ideal Fermi gas which exactly vanishes above pp 
(Fig. [5] inset). This is typical for crystalline solids and 
arises from the Fourier components of the Bloch wave 
functions with G 0. The weights of these components 
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FIG. 5. The momentum distribution of Na determined by a 
DFT-LDA calculation, as compared to the ideal Fermi gas 
step function and to the QMC calculation. Data have been 
rounded to 10""^ corresponding to the accuracy of the calcu- 
lation. 




FIG. 6. Band structure of sodium. Black hue: DFT-LDA 
Kohn-Sham band structure (the experimental lattice param- 
eter at room temperature was used). Red dashed line: dis- 
persion of the ideal Fermi gas evaluated at the corresponding 
density, rs — 3.99. 



(smaller than 0.005 in sodium) quantify the deviations of 
the crystalline wave function from a perfect plane wave 

We can conclude that both band-structure and 
anisotropy effects on the momentum distribution are very 
small in sodium. They affect n{p) by less than 0.02, with 
the maximum value achieved at the level of the discon- 
tinuity. The DFT momentum distribution of sodium re- 
mains close to that of the ideal Fermi gas, significantly 
different from QMC results^ (Fig.El). From Eq. dTO]), we 
obtain the directional averaged Compton profile. Quan- 
titative comparison with the measured Compton profile 
confirms that the bare DFT does not provide accurate 
momentum distributions or Compton profiles (Fig. ([2])). 

IV. CORRELATIONS EFFECTS BY THE GW 
APPROXIMATION 

The Green's function G or the spectral function A en- 
tering into Eq. ([5]) for n(p), can be calculated in the 
Hedin's GW approximation™ where the self-energy is ap- 
proximated by 

I]«^(r,r',..) = * j '^G{vy,u-u')W[v,v\u'), 

(11) 

in terms of the Green's function G and the dynami- 
cally screened interaction W{lu) — e^^{uj)v. Here, v is 
the bare Coulomb potential and e^^{uj) the inverse of 
the dielectric function taken in the random phase ap- 
proximation (RPA).^*^ In an iterative scheme, starting 
with a trial Green's function, G^, for example, obtained 
from the noninteracting band structure or from DFT 
Kohn-Sham orbitals, one calculates the RPA polarizabil- 
ity n^^^ = G^G", followed by the dielectric function 



e = 1 — wll , and the RPA dynamically screened 
interaction M^^^^ = e^pp^{uj)v. With these ingredients, 
the GW self-energy for the first iteration can be cal- 
culated according to Eq. pT|) . and the corresponding 
Green's function is determined by solving Dyson's equa- 
tion, 

G = G° + G°EG. (12) 

In the fully self-consistent GW approximation, further 
iterations up to self-consistency in the self-energy and 
Green's function must be performed. If applied to real- 
istic electronic systems, fully self-consistency is in gen- 
eral a too difficult task, since all functions depend on all 
space-time coordinates. 

Within the so-called GqWo approximation, observables 
are calculated from the self-energy and Green's function 
of the first iteration. The GqWq approximation turned 
oufe^ to be a useful approach for calculating quasiparticle 
electronic structure of realistic systems, yielding band- 
gaps in solids in very good overall agreement with the 
experiment. Different schemes have been explored to ob- 
tain partial self-consistency^2ii^, but fully self-consistent 
solutions are only reported for jellium.^^ Self-consistent 
calculations are appealing since they eliminate the de- 
pendence of the results on the initial trial of the Green's 
function. However, the issue whether a full self-consistent 
GW really improve upon GqWq, remains still controver- 
sial, and likely depends on the observable one is inter- 
ested in. 



A. GW n{p) in jellium as compared to QMC 

The first GqWq electronic structure calculation on the 
jellium model was carried out by Hedin'^^. Lundqvis t^^'^'' 
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FIG. 7. Comparison of our momentum distribution with pre- 
vious results on jellium at Vs — 4. Blue dots: our GoWo 
result; blue solid line: Go Wo result of von Barth and Holm 
(Ref.l4a); red dashed line: GW'^' result of von Barth and Holm 
(Ref . l45fl : black dot-dashed line: GW result of Holm and von 
Barth (Ref.El). 



addressed also the momentum distribution. More re- 
cent calculations, involving self-consistency issues within 
a Gaussian basis set, were done by von Barth and 
Holm."*"*'"*^ Here, we calculated the GqWq momentum dis- 
tribution following the original approach by Hedin J? 

In Fig. [l] we show our calculated momentum distribu- 
tions for jellium for a large range of densities. Start- 
ing from lowest Vs — 1, passing by metallic densities 
Ts =2-5, up to Ts = 10 with gradually increasing cor- 
relation, we observe the departure of the momentum dis- 
tribution from the ideal Fermi gas step function and an 
increasing spill out to higher momenta, associated with 
a reduction of the discontinuity, the quasiparticle renor- 
malization factor Z at the Fermi momentum. We found 
Zpp = 0.86 for Vs — 1, as in Ref. [s^, down to Zp^ — 0.45 
for Ts = 10. 

In Fig. [3 we compare our results to those of Refs. [U 
and li^. Our GqWq n{p) is overall in good agreement 
with the Go Wo result by Ref. 45. The small discrepancy 
at values of p < pi? is not due to the different scheme 
(Gaussian basis set versus sampling of to axis), but rather 
to the fact that for p < pj? we used Eq. ^ instead of 
Eq. ([5]) )^ a point that will be clarified later for the more 
critical case of sodium. 

We then notice that any higher level of self-consistency, 
both on G and W, reduces the level of correlation, and 
thus increases the discontinuity and Zp^ . From compari- 
son with quantum Monte Carlo results of Ref. 12 (Fig.[8|, 
we observe that the GqWq result is already in surprisingly 
good agreement with QMC, especially for metallic den- 
sities. Even for the most correlated density, = 10, the 
QMC discontinuity Z — 0.40 is not far from our GoWo 
prediction of 0.45. On the other hand, any improvement 
towards a fully self-consistent GW has lead to an increase 
of the discontinuity in disagreement with QMC. Notice 



FIG. 8. Comparison of the Go Wo and QMC (Ref. [ij) mo- 
mentum distribution on jellium at several densities. Black 
long-dashed line: rs — 1; red solid line: Vs = 5; dark-green 
dashed line: — 10. 



that Refs. |44| andl45l provide fully self-consistent solutions 
with respect to all degrees of freedom. 



B. GW n{p) in solids 

For Go Wo non self-consistent calculations on realis- 
tic solids, it is convenient^ to start from a Kohn-Sham 



DFT-LDA eigenvalues and eigenfunctions 
zero-order electronic structure 



G"(r,r',c.) = 



0,k(r)</);k(r') 



^ cj - e^k + «77Sgn(eyk - fJ-) ' 



(13) 



used to obtain the screened interaction and the GW self- 
energy, as described above. However, in order to avoid 
double counting of exchange-correlation effects already 
taken into account by the DFT exchange-correlation po- 
tential v^'^, the Dyson equation to calculate the Green's 
function G in first iteration read as 



G = G° + G° (S - v''") G, 



(14) 



In a solid, spectral functions, as well as all other many- 
body quantities, e.g. G and S, can be conveniently de- 
scribed in terms of an orthonormal set of Bloch wave 
functions, 0,yk(r) = 'J2g 't'^k^^''^'^'^^'' ^ already pro- 
vided by the Kohn-Sham eigenfunctions. Expressing the 
spectral function in this basis set. 



A(r,r',w) 



E 



U(r)C'k'(r')^..'(k,k',w), (15) 



the momentum distribution read as 



n(p) 



E 



d{p-k-G) 



duj A^v' (k, k, uj). 

} 

(16) 



o 



In sodium and other nearly free-electron systems, the 
self-energy operator, as well as the spectral function A, 
are almost diagonal on v and i^', and the neglection of 
non-diagonal elements is justified^ in the expression for 
the momentum distribution, 

n(p) = ^5(p-k-G)|0^kl''^.k, (17) 

7iyk = / duj Ai,{k,uj), (18) 

where Aiy(k, w) = Ai/^ (k, k, ui). This approximation does 
not hold in systems where some of the true quasiparti- 
cle wave functions differ from the DFT wave functions, 
and the former expression involving also non-diagonal el- 
ements, A^i^' (k, k, w), should be used. 

In Eqs. (1171) and (|18p . n^]^ provides the correlation con- 
tribution to the momentum distribution from each band. 
For any uncorrelated system, n^k = 1 for all the states 
within the Fermi surface, Ci/k < ^-i^d elsewhere (step 
function) and the strength of the quasiparticle excitation 
remains Z — 1 as for free fermions. However, band- 
structure effects reduce the discontinuity in the momen- 
tum distribution, C, < Z , already for non- interacting elec- 
trons, Eq. Q. Therefore, it is important to distinguish 
Q [the jump in n(p)], from the quasiparticle renormaliza- 
tion factor Z (the jump in rij^k) quantifying the strength 
of excitations at the Fermi surface. Whereas both co- 
incide for a homogeneous system, e.g., jellium, they are 
different in real solids. Within leading order, correla- 
tions induce deviations of rti/k from the step function, 
whereas band-structure effects are already contained in 
the weights inside the summation Eq. (|17l) over all 
bands, reciprocal lattice, and BZ vectors. 



C. GW n(p) by integration on the real a;-axis 

Neglecting off-diagonal contributions in the band in- 
dex, the spectral function, Ec^. (ITSI) . can be expressed 
as 



A^(k,w) 



1 



|3S.(k,L^)| 



TT (c^ - e,k + - 3fiE,(k,w))2 + (3I],(k,w))2 

(19) 

in terms of the band-diagonal self-energy, I]j,(k, cj) — 
Yi^^{k,ijj), calculated within GqWo, the Kohn-Sham 
eigenvalues e^k, and the matrix elements of the DFT- 
LDA exchange-correlation potential 



(20) 



From the spectral function, we can calculate n^k by in- 
tegration along the real axis, Eq. ^TE\\ . Finally, the mo- 
mentum distribution is evaluated from Eq. (jl7p using the 
plane wave coefficients of the underlying DFT-LDA 
calculation. For the GqWo calculation we have used 50 
bands, a cutoff of 5 H (both on the wave functions and 
on the dimension of the polarizability matrices). The di- 
electric function has been calculated on 800 frequencies 
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FIG. 9. Spectral function Av{k,ij) of sodium for the first 
band {v = 1) and k = 0.274 (solid black line); Cumulative 
sum of the spectral function, duj' A{uj') (red dashed line). 
The crossing of the cumulant to the Fermi level provides the 
value of Tifc. The total integral of the spectral function, which 
should be 1, provides the error on n^. A zoom of the A,^{k,ij) 
within the energy interval 3-17 eV is shown by a thin black 
line. 



along the real w axis, up to 2 Ha, and 10 Gauss-Legendre 
knots frequencies along the imaginary axis. With these 
parameters, n(p) at a given p is converged within 10"^. 

In Fig.iniwe show the GqWq spectral function Ay(k, ui) 
for k/kp ~ 0.57. It contains a quasiparticle peak (the 
most intense feature) and some satellites. The distri- 
bution of the spectral weight among the various struc- 
tures can be read off by the integrated spectral function 
(its cumulant). In this case, the quasiparticle peak has 
a weight Z ~ 0.6, low- and high-energy satellites have 
weights of 0.3 and of 0.1, respectively. According to the 
sum rule Eq. ([T]), the total weight of the spectral func- 
tion Ai^CkjUj) is 1. The difference between the integrated 
spectral weight and unity is the numerical error of our 
calculation « 0.01 in this case. The value of rii/k is graph- 
ically provided by the crossing of the cumulant with the 
Fermi energy, i.e., n^i^ = 0.91±0.01 for fc=0.274 and iy=l 
(Fig.®. 

The momentum distribution calculated by the GqWo 
method using the real axis integration is shown in Fig. [TU] 
at the example of sodium. The results are in good agree- 
ment with the QMC calculation over a large range of p. 
When looking at the way n{p) is calculated in a many- 
body approach, Eq. (|18p . the good agreement with QMC 
implies that the GqWo approximation reproduces the cor- 
rect spectral weight repartition between the quasiparticle 
peak and the rest (satellites) . It also implies that the en- 
ergy position of the quasiparticle peak with respect to the 
Fermi energy is correctly reproduced by GqWo- However, 
it does not yet imply a correct energy position of satel- 
lites. Since n{p) is just only sensitive to the satellites' 
spectral weight, GW turns out to be a good approxima- 
tion to describe the momentum distribution, regardless 
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FIG. 10. The GqWo momentum distribution in Na calcu- 
lated along the real (blue squares) and imaginary (violet cir- 
cles) compared to the quantum Monte Carlo (green 
diamonds) data (Ref. Q). In the case of integration along 
the real axis and for k = 0.475 we show both the direct re- 
sult of the integration [n{p) — 0.5] and the corrected value 
of n{p) = 0.78 (the correction procedure is described in the 
text). We show also the momentum distribution in the ideal 
Fermi gas (gray dashed line). 



its ability to correctly describe the satellites' energy po- 
sition. 

Notice that close to the Fermi surface, the straightfor- 
ward integration over real frequencies using a fixed dis- 
cretization grid fails, e.g., the point at n{p = 0.475) = 0.5 
in Fig. llOl obtained from direct integration provides a too 
low value for the momentum occupation. The observed 
underestimation of n^k near the Fermi surface is an ar- 
tifact of the coarse sampling along the w axis. Since 
the quasiparticle lifetime tends to infinity at the Fermi 
surface, the width of the quasi particle peak becomes 
increasingly narrow, and is described by less and less 
points. As we can see from Fig. 1111 the quasiparticle 
peak for k — 0.475 is actually described by only three 
points on the oj-mesh underlying our calculation. As a 
consequence, the integral J dLuA{uj) (black dashed line 
and circles in Fig. 1111) is highly inaccurate. The normal- 
ization of A{uj) turns out to be 0.79 instead of 1, the 
remaining spectral weight, 0.21, is lost due to the under- 
sampling. Using Eq. ^ which involves only the spectral 
function of unoccupied states offers certain improvement. 
However, one should keep in mind that the finite oj-mesh 
employed in calculations on realistic systems might not 
accurately capture all features of the spectral function for 
energies and momenta very close to the Fermi surface. 

The drawback related to the coarse w-mesh can be cor- 
rected by a more appropriated evaluation, without the 
need of recalculation of the spectral function on a more 
expensive fine w-mesh. Based on the assumption that, 
for points k close to the Fermi surface, the (negligible) 
spectral weight of satellites is correctly reproduced by our 
coarse grid, we expand the self-energy inside the spectral 
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FIG. 11. The spectral function for k = 0.475 (close to kp) 
and v = 1. Black continuous line and circles: GqWo spec- 
tral function calculated on a coarse cj mesh; black dashed line 
and circles: integrated spectral weight of the calculated GqWo 
spectral function, /"^ dcj' A{lj'); red dot-dashed line: renor- 
malized integrated spectral weight; blue continous curve: QP 
peak generated by a fit over a Lorentzian; blue dashed line: 
cumulative sum over the QP fit generated curve. 



function, Eq. (jlOp . around their values at the grid-points. 
This justifies the use of a Lorentzian form for the spec- 
tral function around the quasiparticle peak. The spec- 
tral weight is then essentially lost by the poor descrip- 
tion of the quasiparticle peak on the finite lo mesh. This 
procedure corrects the values of momentum distribution 
close to the Fermi surface, in particular, the value at 
n{p = 0.475) — 0.78 is significantly increased compared 
with the value of the direct integration (see Fig. [TU)) . 

In Fig. 1111 we illustrate this procedure where a 
Lorentzian form for the quasiparticle peak is assumed 



1 



QP 



V 



(21) 



and the parameters Zqp, ujq, and ry are determined from 
a fit of the GqWq data such to obtain a total weight of 
1 to satisfy Eq. ([7]). The final result of this procedure is 
the shift of the n{p) point evidenced by the blue arrow 
in Fig. [TTJ The weight above the Fermi level of the fitted 
QP peak (blue line in Fig.[TT]) is reduced, as compared to 
the weight under the the black GqWq points connected 
by direct lines. The obtained integrated spectral weight, 
shown by the blue dashed line in Fig. 1111 provides a value 
of n^fe = 0.80 for = 1 and k = 0.475. The final value of 
n{p) = 0.78 is then obtained as the product of n^k and 
of the wave function coefficient. 



D. GW n{p) by integration on the imaginary tj-axis 

A different strategy to obtain the momentum distri- 
bution consists in deforming the frequency integration in 
the complex plane toward an imaginary w-axis Integra- 
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FIG. 12. Integration contour: the integration along the real 
axis is replaced by the integration along the imaginary axis 
and a quarter circle. 



tion. Therefore, using the normahzation of A Eq. ([7]), we 
write rii,k in the following form 



+00 



+00 



JT-j/k — j4^(k, w) = IH — / (ia;3Gi/(k, a;). 

(22) 

Since the Green's function G,y(k, a;) is analytic in the up- 
per part of the complex plane for uj > fx, the frequency 
integral along the real axis can be now deformed to the 
positive imaginary axis (Fig. [12]). Picking up a contribu- 
tion i7r/2 from the quarter circle at infinite distance, we 
obtain^ 



1 1 f+°° 

2 71" Jo 



(23) 



Thus we can replace the integral along the real cj-axis 
of the imaginary part of G by an integral along the imag- 
inary w-axis of the real part of G, which is in general a 
rather smooth function. The real part of the Green's 
function 5RG can again be evaluated from E solving the 
Dyson equation 



(Ai - e.k + v^^^ ~ + 5E)' ' 

(24) 

In this approach, also the GW self-energy is calculated by 
imaginary frequency integration in Eq. (jlip . The compu- 
tational cost is significantly reduced with respect to the 
integration on the real axis. 

The calculated real and imaginary parts of the self- 
energy T,{uj) along the imaginary frequency w-axis are 
shown in Fig. [13] at the example of sodium. One can see 
that the real part of the self-energy is a much smoother 
quantity on the imaginary axis, in contrast to the imagi- 
nary part^. Close to the Fermi surface, where the imag- 
inary part vanishes, imaginary-axis integration does not 
suffer from the problems encountered in the integration 
along the real axis. So the description of the discon- 
tinuity in n{p) is better in the imaginary- than in the 
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FIG. 13. Real and imaginary parts of the self-energy along 
the u frequency imaginary axis calculated by the non-self- 
consistent Go Wo approach for sodium. We show the real parts 
for three k-points: k = 0.274 (black continuous line); k = 
0.475 (black dashed line); k = 0.6711 (blue dotted line). The 
imaginary part is shown only for k — 0.274 (red continuous 
line and circles). 



real-axis integration (Fig. (TOt . However, undersampling 
artifacts may still occur in the description of states far 
away from the Fermi surface. Here, in contrast to the 
real frequency integration, we have not corrected these 
artifacts, which explains the differences between the two 
integration methods in Fig. 1101 with the real axis inte- 
gration more in agreement with the QMC momentum 
distribution at p far away from the Fermi surface. 



E. Contributions from higher bands induced by 
GW correlations 



In Fig. [14] we analyze the contributions of the different 
bands to the momentum distribution. We show the n^k 
curves, as calculated from Eq. (|23p . for the first v — 1 
band (dotted red lines) and for the second v — 2 band 
(blue double-dot-dashed line), the n^ip) curves includ- 
ing band-structure effects, Eq. p7)) . for the first (red 
dashed line) and second (blue dot-dashed line) bands, 
together with the full n{p) (black solid line). At the non- 
interacting level in sodium only the first band contributes 
to n{p) [see Eq. ([3]), Fig. [5] and [5]. Correlation effects may 
induce contributions from higher bands. From Fig. 1141 we 
observe that in sodium the first = 1 band provides the 
dominant contribution inside the first Brillouin zone, in 
particular for p pf- Correlation effects induce a contri- 
bution n2k from the second band of an almost constant 
0.08. However band-structure effects, taken into account 
by the |0^l^°kP coefficients in Eq. ([T71) . are more intense 
than in the first band (which is almost parabolic and 
planewaves like), and depress this contribution to negli- 
gible values. The contribution of the second band starts 
to be at the level of the first only around p = 1.1 pp 
and for p > 1.3 pp where the li^^^J^^P coefficients are 
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FIG. 14. The momentum distribution n(p) of sodium de- 
termined by a Go Wo calculation with integration along the 
imaginary axis (black squares and solid line). Contribution 
ni (p) from the first (red filled circles and dashed line) and its 
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together with the n2k (blue empty diamonds and dotted line) 
distribution. Lines are guide to the eyes. 



significant. Higher bands {v > 2) provide a negligible 
contribution. 

The discontinuity C in sodium is dominated by the first 
band only. The magnitude of the discontinuity is there- 
fore given by 



(25) 



where Zi,^i_\^p is the renormalization factor of the first 
band quasiparticle peak at the Fermi surface. Since 
l^^i^k P — 0-99j band-structure effects introduce only 
a small 1% reduction of the discontinuity ( in the mo- 
mentum distribution of sodium. 



F. Anisotropies in the GW renormalization factor 

In general in solids the renormalization factor is a func- 
tion of the wave- vector k direction. In Table U we have 
calculated the GqWo Z]^ for the first sodium band, 



^k. = 1 - 



53fiS(kF,cj) 



du) 



(26) 
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along the main symmetry directions in the BZ. In sodium 
anisotropy effects are small also on many-body quanti- 
ties, leading to variations in Z of less than 1 %. The 
GqWq value oi Z = 0.65 ±0.01 in sodium is directly com- 
parable with the Go Wo value oi Z — 0.64 in jellium at the 
sodium average density of ~ 4. Valence electrons of 
sodium look very much like a good realization in nature 
of the ideal homogeneous electron model, both from the 



r-N 
r-p 
r-H 



0.649(1) 
0.643(1) 
0.649(1) 



TABLE I. Anisotropy of the GoWo quasiparticle renormaliza- 
tion factor in sodium. 



point of view of non-interacting, as well as many-body 
interacting observables. We can then compare the GqWo 
value to many different theoretical prediction-^— i^^, in- 
cluding the latest QMC value for jelliumi^, 0.64 ± 0.01, 
which is more accurate than that reported for sodiumM 
(0.70 ± 0.02). The experimental value 0.58 ± 0.07 from 
Ref. [T^ . clearly eliminates several important approxima- 
tions, such as fully self-consistent GW- as well as the 
so-called on-shell RPA approximation. -'^i 



G. GW Compton profile in sodium 

In Fig. [5] we present the GW Compton profile as com- 
pared to the ideal Fermi gas, DFT-LDA, QMC, and the 
IXSS experiment. Both the ideal Fermi gas and the DFT- 
LDA Compton profile appear like an inverted parabola, 
with a discontinuity in the first derivative equal to 1 (or 
nearly for DFT-LDA) at the Fermi momentum. The 
DFT-LDA Compton profile is almost coincident with the 
ideal Fermi gas curve, so that band-structure effects are 
even less appreciable than in the momentum distribution: 
small differences in the two n{p) are further smoothed af- 
ter the integration to get at J{q). 

The discontinuity is clearly reduced in the QMC 
sodium pseudopotential and in GqWq calculations. Both 
QMC and GqWo discontinuities are in agreement with 
the experimental discontinuity within its error bar, 
mostly dictated by the experimental momentum resolu- 
tion and statistical accuracy. The GoWq real axis inte- 
gration Compton profile is practically coincident with the 
Slater- Jastrow QMC result, except at pp where it pro- 
vides a lower discontinuity. A more accurate backfiow 
QMC calculation in sodium would probably further re- 
duce the discontinuity, like it is the case in jellium. As to 
the comparison with the experiment: in principle, most 
possible systematic errors that may influence the exper- 
imental Compton profile, such as the finite experimental 
q-resolution, tend to reduce the experimental value of 
J(0). For this reason, the experimental Compton pro- 
file should be rather regarded as the lowest bottom ex- 
tremum for small q and a highest limiting value for large 



CONCLUSIONS AND OUTLOOK 

We have presented two possible ways within the ab 
initio GqWq approximation to calculate the momentum 
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distribution n{p), the Compton profile J{q) and their dis- 
continuities associated to the quasiparticle renormahza- 
tion factor Zp^. We have analyzed and discussed the 
advantages and drawbacks of both approaches in com- 
parison with QMC calculations and x-ray Compton scat- 
tering measurements on bulk sodium. All the analyzed 
quantities have been found in good agreement with both 
QMC and the experiment. 

In sodium, we have found that n{p) and J{q) are 
mostly determined by the first band and are very weakly 
affected by band-structure and anisotropy effects, in con- 
trast to other alkalines, e.g. Li^"— . This confirms that 
the valence electrons of sodium, in ambient conditions, 
almost perfectly realize the jellium model, even consid- 
ering its electron dynamics^**'^^ . Since jellium is one of 
the most fundamental models to study electronic corre- 
lations, a reedition of experiments on sodium with im- 
proved accuracy may help to clarify open theoretical 
many-body issues. 

Based on the comparison with the Compton profile ob- 
tained by inelastic X-ray scattering spectroscopy exper- 
iments, non-self-consistent GqWq turns out to be a very 
good approximation to calculate the momentum distribu- 
tion and related quantities, while this is not yet evident 
for any attempt including self-consistent GW, both on 
G and W. Therefore, we expect GqWo to provide quan- 
titative reference values for the momentum distribution 
valuable for the development of new exchange-correlation 
approximations within density-matrix functional theory. 
Reverse engineering from GqWq to DMFT, is one of 
the promising ways to improve exchange-correlations in 
DMFT. Similiar strategies have already successfully ap- 



plied regarding improvement of time-dependent density- 
functional theory (TDDFT) by exploiting the Bethe- 
Salpeter equation.—"—. 

Finally, it would be interesting to explore the limits of 
the GW approximation in reproducing the correct quasi- 
particle renormalization factor and the momentum dis- 
tribution as measured in IXSS experiments on other sys- 
tems. An IXSS measure of such quantities, not influenced 
by Hartree or exchange effects like e.g. the bandgap, 
would unambiguously assess the level of correlation in 
these systems and thus the validity of GW in describing 
correlations. A still open question is whether GW can 
describe the quasiparticle renormalization factor in sys- 
tems with increasing level of correlations, up to strongly 
correlated systems. 
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